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Abstract. We compare the energy and count fluxes obtained by integrating over the finite bandwidth of BATSE 
with a measure proportional to the bolometric energy flux, the (p-measure, introduced by Borgonovo & Ryde. We 
do this on a sample of 74 bright, long, and smooth pulses from 55 GRBs. The correction factors show a fairly 
constant behavior over the whole sample, when the signal-to-noise-ratio is high enough. We present the averaged 
spectral bolometric correction for the sample, which can be used to correct flux data. 



INTRODUCTION 

The gamma-ray burst (GRB) spectrum is often peaked 
in the EFe representation (where Fe is the energy flux) 
at some energy ispk in the y-ray band [1]. The limited 
spectral coverage of the used detector might affect the 
assigned measure of the total energy-integrated flux. To 
obtain a true bolometric flux a correction is needed. This 
is most important for the photon flux and less affects the 
energy flux as the energy spectrum is often peaked within 
the detector energy window. 

To circumvent this issue Borgonovo & Ryde [2] pro- 
posed to use the value of EFe at E^^^ as a representation 
of the energy flux, since it is proportional to the bolomet- 
ric flux. We denote this quantity by 9. It was shown that 
the power-law hardness-intensity correlations (HIC), that 
was studied for pulse decays, were better if the (p value 
was used instead of integrating the energy flux over the 
BATSE band. This could indicate that the latter method 
suffers from effects which are not part of the correlation. 
The (p-method thus provides an efficient way of studying 
detailed features in the observed HICs. 

This measure is limited to the cases where the peak 
actually exists within the studied band [^minj^max]- Rep- 
resenting the spectrum by the Band et al. model [1], this 
means that the low and high energy power law indices 
must have a > — 2 and p < —2. In the most common 
case where -E'max 

> (a - p)£'pk/(a + 2), the proportion- 
ality between cp and F=/ FEdE becomes 
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where y^^m = (a + 2)£'min/£'pk, and y^ax = 
(a + 2)£'n-,ax/£'pk- r(a + 2) and P{a + 2,y) are the 
gamma function and the incomplete gamma function, 
respectively (see, e.g.. Press et al. [3]). In the case that 
a and p do not have a strong dependence on fipk, the 
only Epk dependence in ^(a, p,3'min,}'max) is in jmin and 
Jmax- In particular, when the integration is chosen over 
the whole energy range from to 00, there will be no 
dependence at all. The (p-value can be interpreted as 
the integral over the whole energy range or, at least, 
over a range for which y^m < < 1 and y^ax > > 1 , and 
should, under some circumstances, be a better repre- 
sentation of the bolometric flux for the study of the 
hardness-intensity correlation. 



BOLOMETRIC CORRECTION TO THE 
ENERGY FLUX 

For this work, we use a sample of 74 pulses taken from 
55 GRBs, observed by the Burst and Transient Source 
Experiment (BATSE) on the Compton Gamma-Ray Ob- 
servatory. The burst data have a time resolution of multi- 
ples of 64 ms and the count flux is obtained by adding up 
the four energy channels from the Large Area Detectors 
(LADs) of BATSE. We search within bright bursts (with 



peak fluxes larger than 2 photons s 'cm ^) for cases 
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FIGURE 1. Ratio F/(p versus £pk for 846 time-resolved 
spectra from 55 GRBs. The integration of the energy flux F 
was made over the range 20 — 2000 keV. 



containing long, smooth pulse structures with a general 
"fast rise-slow decay". The light curve for each burst in 
the sample was rebinned to achieve a signal-to-noise- 
ratio {S/N) of at least 30. All selected pulses have at 
least 4 time-bins at the given 5 /N. The spectrum for each 
time-bin was fitted with the Band et al. [1] model, allow- 
ing a deconvolution to find the energy spectrum and the 
peak energy. The energy spectrum was integrated over 
the energy band of the BATSE detector (20 - 2000 keV), 
for the strongest illuminated LAD, to find the instanta- 
neous energy flux. 

In Figure 1, the spectral bolometric correction for 846 
time-resolved spectra from the 55 GRBs is shown. Av- 
erage values can be used to correct the energy flux data 
when = £'pk(0 is known. The data show a greater 
dispersion for smaller E^^. This is most Ukely a bias in- 
troduced by the fact that most of the spectra correspond 
to the long decay phase of pulses, during which the inten- 
sity is positively correlated with the hardness. Therefore, 
the signal-to-noise ratio is usually much higher at the 
peak of the pulses when the spectra are harder, i.e., for 
higher E^^, and decreases to the chosen Umit {S/N = 30) 
as the pulse decays. 



BOLOMETRIC CORRECTION TO THE 
COUNT OBSERVATIONS 

We now redo the above analysis on the count flux 
instead. The count data have not been deconvolved and 
for proper physical interpretation the effective correction 
must be estimated and understood. 

Figure 2 shows the ratio of the instantaneous energy 
flux F (in keV s~^cm~^) to the corresponding count flux 
Fc (in counts s~^) for the same time-resolved spectra as 
in Figure 1. The data were fitted with a second-order 
polynomial which was found to be 0.05235 -|- 1.924 x 
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FIGURE 2. Ratio of the instantaneous energy flux to the 
count flux versus the corresponding E^\^ (in keV), for the 
same sample of time-resolved spectra as in Fig. 1. A best- 
fit second-order polynomial 0.05235 + 1.924 x lO^'^^pk - 
6.624 X 10~^£'pjj^ is also shown. The dispersion parameter is 
0.1516. 
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FIGURE 3. Ratio of the (p-measure to the count flux Fc ver- 
sus Spk. A second-order polynomial fit to the data 0.01597 4- 
8.320 X 10"5£pk - 2.829 x W'^E^^ is also shown. The disper- 
sion parameter is 0.07605. 



lO-^^Epk - 6.624 x lO^^E^y, (with £pk given in keV). The 
dispersion was measured as the mean of the squared ra- 
tios between the fit residuals and the fit expected values. 
The F /Fc fit has a dispersion of 0.15. 

In Figure 3, the ratio of the (p-measure (in keV 
s^'cm^-^) to the count flux (in counts s^') over the same 
sample of spectra is shown. Again the data were fit- 
ted with a second-order polynomial given by 0.01597 + 
8.320 X IQ-^Ep^ - 2.829 x 10~^E^^ (wifli % in keV). 
This adjusts the data with higher accuracy for small ispk 
values. The calculated dispersion is 0.08, significantly 
smaller than in the F/Fc case, shown in Figure 2. 



DISCUSSION 



To interpret the BATSE spectral data, one must have in 
mind that due to the GRB spectra being broad compared 
to the observed energy window, one must apply a bolo- 
metric correction factor which is unknown. To circum- 
vent this problem we propose the use of the cp measure 
that is also proportional to the bolometric flux. 

The flux F is found from integrating the deconvolved 
spectrum. This spectrum is model-dependent as the de- 
convolution is based on the model spectrum. This could 
introduce additional scatter into the determination. A 
second problem, when one aims at studying single pulses 
in the light curve, arises from the fact that the observed 
spectra may contain contributions from other pulses 
which even could be unresolved. Furthermore, additional 
soft components [4] could also affect the measured flux 
value and thus weakening the correlations. This makes 
the (p-method better to use. 

The F/(p ratio shows an approximately constant be- 
havior over the studied sample, with a larger dispersion 
for smaller Ep^. This fact can be explained as a bias in- 
troduced by the HIC correlation. 

We also analyze the correction for the count flux. 
The (p/Fc ratio is shown to have a substantially smaller 
dispersion than the F/F^^ case. Therefore, the effective 
correction for the count flux may be better estimated. 
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